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Abstract 

We study relaxation processes in spin systems near criticality after a quench 
from a high-temperature initial state. Special attention is paid to the stage 
where universal behavior, with increasing order parameter m(t) ~ t e , emerges 
from an early non-universal period. We compare various algorithms, lattice 
types, and updating schemes and find in each case the same universal behavior 
at macroscopic times, despite of surprising differences during the early non- 
universal stages. 
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Temperature quenches in spin systems and the ensuing relaxational processes have been 
a much-studied subject in the recent past [0-0]. Especially for a quench from a high- 
temperature initial state to the critical region, a new universal regime was predicted by 
Janssen et al. ||. In its most pronounced form this phenomenon, termed universal short- 
time behavior (USTB) in the following, occurs in a model with purely relaxational dynamics 



model A according to Ref. f20|). Starting from an initial state with T>T C and a small 



initial magnetization mo, the magnetization increases as 

m(t)~m t e (1) 

for a macroscopic time span before it reaches a maximum and eventually decays to the 
equilibrium value zero. 9 is an independent new exponent determined by the non-equilibrium 
initial state that cannot be expressed in terms of equilibrium exponents. Its value in the 
two-dimensional Ising system for instance is 9 ~ 0.19 P,|T5|. 



The results of Janssen et al. 0|| allowed the interpretation of earlier Monte Carlo (MC) 
simulations Q], and later the power law ([!]) was also directly verified in the "computer 
experiment" . Further it was pointed out by Li et al. || that the USTB can be exploited to 
determine equilibrium exponents from the early non-equilibrium. In the sequel, this method 
was further developed and applied to a number of systems [0|§. Moreover, USTB in other 
dynamic universality classes |§, in systems with a tricritical point |l0j, and in disordered 



and dilute spin systems was studied [11]. Finite size effects were analyzed in Refs. [12,13 



and the USTB near surfaces was studied in Ref. ||14|| . Further, the close relationship between 
USTB and "damage spreading" was pointed out flI5| , USTB in a different context, for quasi- 
long-range order evolving after a quench to the Kosterlitz-Thouless phase, was investigated 
TEL and a general scaling invariance in the short-time regime was found [|17]]. Possibly 



also related to USTB is the "overshooting" of the order parameter beyond the equilibrium 
domain magnetization for quenches below T c |Tj|], and the issue whether there exist still 



more independent exponents for relaxational processes was raised in a recent preprint JT9[ 

A simple physical argument for the growth of the magnetization in ([!]) was given by 
Janssen |||21|]. Consider a system that is quenched to some final temperature Tf (not 
necessarily T c ), again with initial magnetization m®. Then for Tf <C T c , m(t) should grow 
after the quench, towards the equilibrium value selected by itlq. If in contrast Tf S> T c , m(t) 
is expected to decay to zero rapidly. Hence, there should be a limiting temperature T z where 
the qualitative behavior changes. 

As the initial correlations are short-ranged, the natural candidate for 7] is the critical 
temperature of the mean-field (MF) theory T C MF , and with the real T c of spin systems being 
always smaller than T^ F , it would be an immediate consequence that m(t) increases for a 
quench to the critical point. However, as argued in Ref. it is not possible to derive 



the power law (|1|) from this scenario. The power-law growth is rather a phenomenon that 
occurs when the time-dependent (growing) correlation length £(t) has become macroscopic, 
i.e., much larger than the lattice spacing a (compare Fig. 1 below). The derivation of (|TJ) 
is thus beyond the scope of MF theory. 

So far numerical investigations have been mostly carried out with the heat-bath (HB) 
algorithm p2|j23|] . A comparison between HB and Metropolis (ME) algorithm |22| was 



performed for the Potts model by Okano et al. ||, and it turned out that concerning the 
universal behavior both algorithms yield compatible results, but differences occur for early 
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times. 

The main purpose of this Letter is a more systematic examination of the issue of univer- 
sality. Are equilibrium exponents determined with the USTB really independent of factors 
like the algorithm (HB or ME), the updating scheme (random or sequential), and the lattice 
type (nearest or next-nearest neighbor coupling, square or triangular lattice)? How does 
universal behavior in the regime with £(t) 3> a emerge from the non- universal early stage 
with ~ a? And closely related: Is it really a mean- field ordering process during the 
microscopically early stages, or has this simple picture to be refined? 

We answered these questions by solving the master equation for early times (during the 
first single-spin updates) as well as by MC simulation for later times. Our work reveals a 
number of interesting and surprising details about algorithms for critical dynamics and puts 
the USTB as a method to access equilibrium properties from a non-equilibrium (though 
universal) regime on a much firmer basis. 

Information about thermal averages in the early non-universal stage of the relaxational 
process can be obtained by solving the master equation 



^^ = J2[W( S '^ S )P( S ',t)-W(s^ S ')P(s,t)] , (2) 

s' 

where s denotes a spin configuration, W(s' — > s) is the transition probability, P(s,t) the 
probability to find configuration s at time t, and the sum extends over all possible con- 
figurations. The analytic integration of Eq. (0) for large systems of coupled spins is not 
feasible. However, for the high-temperature initial state consisting of N uncorrelated spins 
in a magnetic field H, characterized by the Boltzmann factor 

N 

P(s, 0) = Z X exp(# s i) with Z o = ( 2 cosh H ) N , ( 3 ) 

i=l 

the analytic treatment for very early times, £ <C 1, is possible. (The time is expressed in 
units of MC steps per site (MCS).) 

Consider the first single-spin update after the quench. It takes place in an environment 
that is coupled to a heat bath at the final temperature Tf, and H is switched off. Decisive 
for the very early stage is whether after the first update on average the magnetization is 
reduced or increased. The respective tendency survives as long as the system still closely 
resembles the initial state, i.e., as long as the number of single-spin updates is much smaller 
than N. 

Without loss of generality one may choose spin 1 to be updated. From (§) one straight- 
forwardly derives Am := m(£ = 1/N) — m to be 

Am = 2 N- 1 [W( > +) P(- 0) - W(+ -> -) P(+, 0)] , (4) 

s 

where the W's are the probabilities for spin 1 changing sign with all other spins remaining 
unchanged. Especially for small H (corresponding to small mo), we find from (§) the simple 
result 

Am = -4HN- 1 Y / W{+^-) 11+ £ sA , (5) 
s V ielN ) 
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where the second sum extends over the interacting neighbors of spin 1. 

From (||) we calculated Am and the limiting temperature T/ for the Ising model on a 
square lattice with nearest-neighbor interaction J for HB and ME algorithm. In d = 2 
(the system that will be also studied by MC simulations below) the explicit results are 
T t HB = 3.0898 ... and T t ME = 1.5885 ... for HB and ME algorithm, respectively. (Temper- 
atures are expressed in units of J/ks-) For comparison, the (exact) critical temperature is 
T c = 2.2691 . . . and T C MF = 4. Hence, for the HB algorithm T\ is indeed above T c , while, 
surprisingly, with ME even at T c the magnetization decays t < 1. 

We calculated T; also for other dimensionalities. In the limit d — > oo (where the number 
of nearest neighbors becomes infinite) one finds TJ —>■ T^ F for both algorithms considered 
above. For d = 1 both yield 7] = 0. 

From this analysis we conclude: First, the limiting temperature in general depends on 
the algorithm. Only for d — > oo it turns out to be T C MF . Second, we are left with the 
puzzling result that for the ME algorithm the limiting temperature lies even below T c , and 
therefore one would not expect to see an increase of the magnetization at T c . In any event, 
the simple explanation that the non-universal stage preceding the USTB is a MF ordering 
process is in general not correct. 

In order to learn about later stages, t ^ 1, especially the crossover from microscopic to 
macroscopic behavior, we had to resort to MC simulations. These were carried out for an 
Ising system on a square lattice with linear dimension L and periodically coupled boundaries. 
Single spins were randomly selected and updated. In order to obtain thermal expectation 
values we generated a large number of histories, each starting from a new initial configura- 



tion, and calculated mean values [E2 




Fig. 1: Three snapshots of the temporal evolution of a spin configuration for L = 300 and m,Q = 0, 
generated with the HB algorithm and random updating. Displayed in each picture is half of the 
system. From the visual appearance there is no difference HB and ME algorithm. 
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Snapshots of the temporal evolution of a single configuration for a square lattice with 
N = 90 000 spins are displayed in Fig. 1. The left picture shows the initial state. Next to 
it the configuration after 2 N updates corresponding to t = 2 is depicted. At this point the 
average domain size and (with that) the correlation length are already substantially larger 
than the lattice spacing. This is the stage where universal (macroscopic) behavior emerges 
from the non-universal (microscopic) regime as discussed in more detail below. At t — 100 
the correlation length is of the order of the lattice size L. 

Results for m(t) at Tf = T c with L = 20 and tuq = 0.05 are displayed in Fig. 2. The 
HB curve (solid line) monotonously increases and is consistent with a power-law for t > 1.5. 
In the case of the ME algorithm (dashed line) the behavior is qualitatively different. As 
expected from the our analytic results, m(t) indeed drops initially, but has a minimum at 
t ~ 0.3, and then increases to assume the power-law form for t ^ 2. Thus, despite of the 
anomalous time dependence of the ME curve in the non-universal regime, for macroscopic 
times it agrees with ([I]). This is in accord with the findings of Okano et al., where sequential 
updating was used and, thus, the details of the temporal evolution were not uncovered. 
Later the profiles in Fig. 2 have a maximum and then decay to the equilibrium value zero 
23. 
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Fig. 2: Order-parameter profiles for L = 20 and ttiq = 0.05 obtained with HB (solid line) and ME 
algorithm (dashed line) in double-logarithmic representation. The small diagram inserted shows 
the data for small times in double-linear representation. 

Taking into account these results, the natural question to ask is whether there exists 
an algorithm-independent limiting temperature, a "dynamic MF temperature", where the 
macroscopic behavior changes from increasing (USTB) to decaying. It turns out that for the 
HB algorithm this limiting value coincides with T^ B \ the profiles have mostly one extremum. 
In order to determine the corresponding limit for ME, we generated a number of profiles 
for temperatures above T c , seeking the one that shows a saddle point. We detemined the 
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corresponding temperature as 2.70 (2). This number did neither depend significantly on the 
system size nor on m . However, it does not agree with the corresponding value of the HB 
dynamics T l HB ~ 3.1. 

Eventually we compared magnetization profiles in the universal regime for a system 
with L = 40 and mo = 0.03 at the critical point, Tf = T c , for various combinations of 
algorithms, lattice types, and updating schemes. For the square lattice with nearest-neighbor 
interactions we combined ME and HB algorithm with random and sequential updating (4 
curves) , for the triangular lattice with nearest-neighbor interactions we used both algorithms 
and sequential updating (2 curves), and for a square lattice with additional next-nearest- 
neighbor interactions we used the HB algorithm and sequential updating (1 curve) [25|| . 

As can be seen already from Fig. 2, even though the initial power law is assumed by both 
profiles depicted there, the heights and locations of the maxima depend on the details of the 
method, besides the differences for early times. However, in all cases studied it turned out 
to be possible to map the data onto a single curve for times t > 20, by constant rescalings 
both axes. The result is shown in Fig. 3. On the semi-logarithmic plot (small insertion) 
the individual profiles cannot be distinguished. When both axes are plotted logarithmically, 
on the other hand, the short-time regime is more pronounced, and significant differences 
for t < 10 become visible. The pure power law (solid line above the data) is plotted for 
comparison. 
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Fig. 3: Data collapse of seven magnetization profiles obtained for L = 40 and m,Q = 0.03 for dif- 
ferent combinations of lattice types, algorithms, and updating schemes (as described in the text) in 
double-logarithmic representation. The ME results for the square lattice with sequential (dashed 
line) and random (dotted line), as well as for the triangular lattice with sequential update (dashed- 
dotted line) are singled out. The small diagram inserted shows the same data in semi-logarithmic 
representation. The power law (Q) with 9 = 0.19 is plotted for comparison. 
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Singled out are the ME curves, for the square lattice with sequential (dashed) and random 
(dotted) updating, and for the triangular lattice with sequential updating (dashed-dotted). 
In particular when the ME algorithm is combined with sequential updating, the power law 
is assumed only for t ^ 10. This is consistent with the findings of Okano et al. || for the 
Potts model, and can now be interpreted as a consequence of the anomalous behavior of the 
ME algorithm for t ^ 1. The HB results assume the power- law form much earlier. Most 
importantly, however, all seven curves show the USTB as expressed in ([I]) for later times, 
*>10. 

In conclusion, we studied the short-time behavior in relaxational processes after a tem- 
perature quench. We found surprisingly different temporal evolutions with heat-bath and 
Metropolis algorithm during early, non-universal stages. Nevertheless, the characteristic 
short-time power law ([!]) turned out to be a rather robust phenomenon, occurring indepen- 
dently of the algorithm, the lattice type, and the updating scheme, provided the systems be- 
long to the same dynamic universality class — the Ising system with short-range interactions 
and non-conserved order parameter was the example studied above — and the correlation 
length has grown substantially larger than microscopic scales. 
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